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ABSTRACT 


The pressure drop characteristics of centrifugal pumps during pressurized startup 
were studied analytically. Mathematical expressions for negative head characteristics 
were developed and applied to the prediction of pressure drop characteristics for three 
pumps. Comparison of calculated and experimental performance points verified the 
analytical method developed. 
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SUMMARY 

The pressure drop characteristics of centrifugal pumps during pressurized startup 
were studied analytically. A method for predicting pump pressure drop characteristics 
was developed. Values were calculated for three pump configurations and compared 
with experimental data. Enough test data were available for the two pumps with vaned 
diffusers to verify the analytical method quantitatively. In the third case, a pump with a 
vaneless diffuser, the procedure was verified qualitatively. 

The method developed is based on the analysis both of work done on the pump’ s im- 
peller and of the hydraulic losses within a pump when it is subjected to flow forced by a 
pressurized tank. An empirical correction factor based on an analytical extrapolation 
of the data from the pump's high-speed characteristic is used to adjust the analytical 
equations for the generally unknown losses due to the mismatch of flow and blade angles 
at the entrance of the diffuser. 

The close agreement of calculated and experimental performance points verifies the 
analytical method. The study demonstrates that it is possible to predict the pressure 
drop characteristic of centrifugal pumps with vaned diffusers provided that the high- 
speed characteristic is known. To predict the pressure drop characteristic of a centri- 
fugal pump with a vane less diffuser, knowledge of the high-speed characteristic is not 
required. 


INTRODUCTION 


In the development of the chemical and nuclear rocket engine systems, it is neces- 
sary to know system performance not only at the design point but also in a wide range of 



the off -design conditions. The need for this information is particularly severe in the case 
of engine systems where the hydraulic losses in the pump during startup affect the engine 
transient characteristics. 

In the case of the axial turbopumps, methods (ref. 1) are available to calculate pump 
characteristics in a range 80 to 120 percent of the design flow point. Also, an analytical 
method has been developed, verified by experiments, and published (ref. 2), for the pre- 
diction of pressurized startup performance of a multistage axial pump. 

Prediction of the off -design operation of centrifugal pumps presents a more difficult 
problem, especially for pumps equipped with vaned diffusers. The difficulty lies in the 
lack of information on hydraulic losses at the discharge of centrifugal impellers and at 
the entrance to the diffusers. 

The usual method of evaluating the performance of a centrifugal pump at and near the 
design point is based on a comparison of the design parameters with the parameters of a 
similar pump of known performance. For strong off -design conditions such as those near 
the shutoff or maximum capacity points, it is necessary to correlate the empirical data 
from tests on a diffuser similar to the diffuser used. For extreme off -design operating 
conditions encountered during pressurized startup, there are no diffuser loss test data 
and no performance prediction method is available. 

In this report, analytical expressions are developed to describe pump performance 
during pressurized startup. These analytical methods are applied to the following pump 
designs: (1) a radial centrifugal pump with a vaned diffuser, (2) a centrifugal pump with 
an impeller blade exit angle of 60° and equipped with a vaned diffuser, and (3) a centrifu- 
gal pump with an impeller blade exit angle of 35° and equipped with a vaneless diffuser. 

The pressure drop characteristics are calculated and compared with data from tests 
run on the first two pumps. An additional qualitative verification of the theory is obtained 
using the design data of the third pump. 


APPROACH 

The generalized characteristic of a pump extends from positive to negative pressure 
differences. A study of complete generalized characteristics is given in reference 3. 
Figure 1, derived from this study, shows the more important features of the generalized 
characteristics of turbopumps. The generalized variables are AP/N^ and w/N. (All 
symbols are defined in appendix A. ) There are three regions of pump operation during 
pressurized startup. Region A is called "the reversed turbine operation” region; the 
pump operates as a hydraulic turbine with the torque supplied by fluid fed from a pres- 
surized tank. Impellers with a blade discharge angle other than 90° will act as a hydrau- 
lic turbine when subject to induced tank flow. The characteristic of such a hydraulic 
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Figure 1. - Generalized pump characteristic extending into posi- 
tive and negative pressure difference. 


turbine can be calculated or can be obtained by tests. Region B is called ’’the energy 
dissipation region” where the torque is supplied by the flowing fluid and the primemover 
simultaneously. In this region, the pressure difference across the pump is still nega- 
tive. Because the degree to which the primemover supplements the fluid torque is 
unknown, the characteristic in this region cannot be calculated. Region C is the pump 
operation region where a torque is applied to the pump by a primemover and there is a 
positive pressure difference across the pump. In this region, the pump characteristic 
may be approximately calculated or established by experiments. 

There are two points in this figure which must be discussed. Point 1 is the pump 
zero -efficiency point. A pump reaches this point when the applied torque serves to 
balance the losses within the pump and no positive pressure difference is generated. 

Point 2 represents turbine (reverse hydraulic turbine) zero -efficiency point, that is, the 
zero-torque point or windmilling. At that point the impeller rotates at a speed that 
corresponds to the available pressure drop (inlet and outlet conditions) and is consistent 
with the flow that can be handled. At this point, no useful accelerating torque is availa- 
ble. The centrifugal pump characteristic (region C), to the right of the design point, has 
an increasingly negative slope. The hydraulic turbine characteristic (region A) is much 
steeper, and its slope varies little with the flow parameter w/N. The two characteris- 
tics merge in region B where an adjustment of slopes must take place. 
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Figure 2. - Pump performance during pressurized startup for 
analytical model. 


A more convenient form of the generalized characteristic for the startup calculation 

. o 

is shown in figure 2. The new variables are AP/w and N/w. The known pump char- 
acteristic and the calculated turbine characteristic are drawn as solid lines. The dashed 
lines represent analytical extrapolation of the pump characteristic beyond its experimen- 
tal range and into regions B and A, and extrapolation of the turbine characteristic into 
region B. It is required that the turbine characteristic meet the pump characteristic at 
point 1. Because extrapolations are used, the two lines may not actually meet and, in 
practice, the requirement is that the two extrapolated characteristics must approach one 
another with the minimum deviation that can be found. 

The path of the startup process between points 2 and 1 is not shown because it is es- 
sentially empirical; it should lie between the two dashed lines. Point 3 represents the 
locked rotor pressure drop or the pressure drop corresponding to the very initial mo- 
ments of a startup process before any rotation is obtained. The location of point 2 de- 
pends on pump design. For an axial pump, point 2 is very close to point 3, as shown in 
appendix B where, from reference 2, the key results of a startup performance analysis 
for an axial pump are given. For a radial centrifugal pump, point 2 will be at point 3. 
For other designs, it will be anywhere between points 3 and 1. The ability to calculate 
the turbine characteristic or the zero -speed point for the radial pump makes it possible 
to define the complete pump startup characteristic. 
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DEVELOPMENT OF THE ANALYTICAL METHOD 


Data Requirements 

To calculate centrifugal pump startup characteristics, the following information must 
be available: 

(1) The pump geometry including blade angles 

(2) The design point flow, speed, and pressure rise 

(3) The pump characteristic in a range within approximately ±10 percent of the de- 

sign point flow. 

To calculate the startup characteristic of an axial pump, the pump characteristic was 
not required; the velocity triangles at the design point were sufficient (ref. 2). 

The difference between the analytical approach used in calculating the pressure loss 
characteristic of the axial pump of reference 2 and the approach needed for a centrifugal 
pump lies in the inability to predict losses resulting from the mismatch of the flow and 
diffuser blade angles of a centrifugal pump. However, with a known pump characteristic, 
the calculated startup characteristic may be adjusted for those losses by the use of a 
correction factor. 


Head Characteristic 

The energy equation for an incompressible and nonconducting fluid and for a steady 
adiabatic turbulent flow, where the change in potential energy due to change in height is 
neglected, is (ref. 4) 


Ah =^+^W + ^cp (1) 

2g g g 

which states that the change of the static head equals the change of the kinetic energy 
(based on the average conduit velocity), the work done on or by the impeller W, and the 
viscous losses cp (the energy that is dissipated due to viscous effects per pound of 
fluid). Equation (1) expressed in terms of the total head change becomes 

Ah = — W n — — cp (2) 

g g 
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The low -velocity hydraulic turbine flow criteria and a one -dimensional analysis, sim- 
ilar to the approach used in reference 5, are used to describe the flow dynamics of the 
system under consideration. Consequently, leakage, windage, and mechanical losses are 
neglected and the work done by the fluid on the impeller is considered negative. 

To expand equation (1), use will be made of the fluid velocity components over blades. 
For simplicity, in subsequent equations, the parameter g c will be replaced by g. Fig- 
ure 3 gives station designations used in pump analysis. 



Change in Kinetic Energy 

The change in kinetic energy is completely described by the pump inlet and the pump 
outlet velocity head difference -(v^ - V?^2g where V e is the velocity at the pump dis- 
charge f lange and V. is the velocity in the inlet pipe. 


Work 

The work done by the fluid on the impeller must be described by the blade character- 
istics. The work done on the blades, by Euler's equation, is -(u^V u ^ - UgA^ y)/z. 
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Frictional Losses 


The losses resulting from flow through the stage are of two types, losses due to fric- 
tion and losses due to turbulence induced by the mismatch of blade and flow angles caused 
by the off -design operation. The frictional losses will be calculated using impeller and 
diffuser exit velocities. Actually, average fluid velocities in the flow passages should be 
used. Because the frictional losses during startup (low flows) are small compared with 
all other losses, the resulting error will be negligible. The frictional losses for the im- 
peller and for the diffuser, respectively, are 





-4f 



d h 2 g 




2g 


where c is the mean-diameter blade length and d^ is the hydraulic diameter of the 
fluid passage formed by the average blade spacing s and the average blade height b. 
The friction coefficient f will be taken equal to 0. 004, an average value over a wide 
range of Reynolds numbers (ref. 4), because as noted, the frictional losses are small. 


Turbulence Losses 

The turbulence losses which occur at the inlet of the impeller and the diffuser may be 
derived by analyzing the impeller inlet, impeller outlet, and diffuser inlet velocity 
triangles. 



Figured - Impeller inlet velxity 
triangle at off-design operation. 
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An efficient operation implies that the fluid angle a\ will equal the design angle 

t t » A 

c*! ( j. At off -design conditions, + 01 ^ Consider an arbitrary impeller inlet veloc- 
ity’ triangle XYZ as shown in figure 4. The fluid approaches the blades with a relative 

t t t t 

velocity at an angle ot^. Because ^ a tangential velocity component of 

adjustment ZZ' will be required to force the fluid to follow the blades. The energy loss 
resulting from this change of direction equals -(ZZ’) 2 /2g. From the geometry of the 
trapezoid XYZ’Z 


sin(ci! 1 . + a.) 

ZZ’ = u i- V x = u x - k r V 1 


sm a 


l,d 


where 


k 


r 


sin(o! j d + 


sina l,d 


and the loss is 


2g 

Similarly the diffuser entrance loss may be evaluated in terms of the impeller outlet re- 
lative velocity V 2 and of the tip blade speed u 2 , corrected for the diffusion effects in 
the annular space between impeller tip and the diffuser base circle. The velocity trian- 
gle at station 2 will be recalculated for station 3 by multiplying the tangential velocity 
components and the axial velocity by ratio r 2 /r 3< 

The diffuser entrance loss will be 


(-3 - k s V 3> 2 
2g 

which is -(u 0 - k v’ 9 ) 2 /2g adjusted for free vortex diffusion between stations 2 and 3. 
The loss coefficient is 
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_ sin (“3 + a! 3,d ) 

K s . 

sm a 3jd 

By a simple trigonometric transformation it may be shown that k r = u^/V and 
k = Ug/Vg which indicates that turbulence loss calculation assumes no loss at the pump 
design point operation and assumes a definite incidence angle at which this occurs. While 
this is of no consequence with the impeller inlet loss, it is generally not the case with 
the diffuser inlet loss. 

Experimental data reported in appendix C indicate that diffuser inlet losses may be 
important. In this experimental work, generally the diffuser loss is never zero and the 
minimum loss, depending on design, may be significant. Furthermore, the incidence at 
which this minimum occurs depends on pump and diffuser blade configuration. Examina- 
tion of the design parameters of the diffusers under consideration indicate that the inci- 
dence related to the minimum loss and the corresponding pressure drop cannot be corre- 
lated with the design data. A successful correlation would require an extensive study of 
different designs. For the same reason, the losses at off -design conditions are even 
more difficult to evaluate. A further complication arises from the fact that pump design 
point (best efficiency) seldom coincides with the minimum diffuser inlet loss. 

This lack of information on the diffuser entrance loss requires that an empirical 
correction factor X must be used to adjust k g to a new loss coefficient K g = Xk g . 

Stage Exit Loss 

It will be assumed that except for vf/2g all velocity head at discharge from the 
diffuser (vaned or vaneless) is lost. The loss will be equal to (V^ - V e y 2g. 

Generalized Negative Head Characteristic 

By use of the relations obtained for the work done and for the various losses, 
equation (1) becomes 

2g Ah s = V? - v| - ly\ - vl) - [(kj ; S V l ♦ k, r vj 2 ) + 2(u lVu> , - u 2 V u _ 2 ) 

+ (», - k r V,) 2 + <u 3 - k s v’ 3 ) 2 ] (3) 
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To express the head loss in terms of and only, the following relations, de- 
rived from the velocity triangles will be used: 

Vj cos 

.0 for = 90° 

V u,2 = u 2 -Vgcos 

O O I o t t 

Vg = Ug + V 2 - 2 u 2 V 2 cos Q ! 2 
V u, 2 r 2 = V u, 3 r 3 

t t 

In addition, it will be convenient to eliminate V 2 by writing V 2 = yV^ where y is 
a proportionality constant. The two velocities V 2 and are always proportional to 
each other because for any given slip their orientation is fixed. Now since N = 30u/7rr 
(in SI units, N = u/r) and by continuity equation Q = sin a^, the expression for 

negative head characteristic may be written in measurable quantities such as N and Q 
instead of the corresponding velocities u^ and V^. 

After substituting and rearranging, the head parameter Ah g /Q^ of the pump is 
expressed as a function of the speed parameter N/Q and a number of design parameters 
as follows: 



Ah 


f, s 


2 , 

y k . 

J f , r 


Q 2 2g^ 2 2g^ 2 1 8g7r 2 r^b^ sin 2 a> 4 2g ^ 


2222 22 22 
8tt- gr^b^ sin a 4 2g^ 1 sin a ^ 2gr^ sin 


2 , 2 / 2 2 * 
r 2\ k s A C ° S + 1_ 


ITT. 


sin 2 o'! 47r 2 r^b 2 / 30 S«1 sin «1 


cos a ^ + y 


cos a 2 - kp - sin ( — 


V 2\( V 2 


x k„ 


, 2 2 T 
y cos o ?2 i 
■ + 


ft 2 sin 2 4jrr 2 b 2 


2 2 

N\ . 71 r l 


Q/ + 1800 g 


- 2 /— 


r l/\ r 3 


(4a) 
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r 


In special cases, this equation simplifies greatly. If the diffuser is located in the imme-' 
diate vicinity of the rotor tip rg = r^ and no diffusion need be considered; for a pump 
equipped with a vaneless diffuser k g = 0 and also r^ = rg. In the SI system of units, 
this equation becomes 


1 f _1 j_\ J k f,s , ? k f,r , k r 

Q 2 2 gft 2 2 g£ 2 2 \ 8 g 7 r 2 r^b 2 sin 2 a 4 2 8 g 7 r 2 r^b 2 sin 2 a 4 2 gft 2 sin 2 2g£2 2 sin 2 


'r 0 \“ k* /y^ cos 2 oi< 


2g sin 2 a 1 


47 r 2 r 2 b 2 / g^i sm «i 


/ r 2\ 

cos q ? 1 + y / — j cos a 2 ” k r “ 


^2 ’ K r " s ^ n i — 


x k„ 


2 2 
fy cos o! 2 


9.2 „ 2 2,2 

sin 477 1*2^3 


2 

N\ + ^1 


Q/ 2g 



(4b) 


Equations (4) are applicable to single-stage centrifugal pumps. If a multistage cen- 
trifugal pump is considered, all terms in these equations, except the first two, will be 
written for each stage using appropriate velocity triangles. 

Equation (4a) is the final general expression for the negative head characteristic of 

centrifugal pumps in terms of head difference Ah and volumetric flow Q. For practi- 

s 

cal applications, it will be more convenient to use pressure drop and speed parameters, 

• 2 

AP/w and N/w, respectively. 

Equatiom(4a) will be applied directly by assuming that there is no slip. This assump- 
tion may be justified as follows: in radial centrifugal pumps the fluid will discharge from 
the impeller at the blade angle of 90° with no tangential component capable of producing 
a starting torque. For a blade discharge angle other than 90° the initial torque will tend 
to rotate the pump but will be opposed by friction. The mechanical losses such as bear- 
ing friction, not considered in the analysis, may however, prevent the rotation. Once 
the rotation is set in, the slip, if any, will be small and may be disregarded. 

The blockage due to the number of blades must be accounted for. It will influence 
the shape and size of the velocity triangle at the impeller discharge and must be carefully 
calculated. 


Point of Zero Work Done on Impeller 

By equating to zero the work term of the negative head characteristic and by making 
proper substitutions and rearrangements, the value of the speed parameter N/Q for the 
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zero -work point or the zero -torque point as it will be called, may be calculated from 


(ujV . - u 9 V 9 ) = ^ (r 1 cos a. + r„y cos aJ 

X U, X £t U, u sin Q!j 


N = 0 


30 


(5a) 


which, for axial inlet to the rotor (a^ = 90°), simplifies to 


\t 30g cos a 0 

N = 2 (5b) 

Q 7712^2 

In the SI system of units these two equations are unchanged except that the factor 30A 
drops out because it is contained in N. For radial blading when slip is not considered, 
a 2 = 90° and N/Q = 0; that is, there will be no windmilling and a torque from the tur- 
bine is necessary to start the rotation. 


VERIFICATION OF THE ANALYTICAL METHOD 


Mark III - Model IV Turbopump 


Data. - Based on experimental data, the high-speed pump map (unpublished data from 
Aerojet -General Corp. ) is shown in figure 5(a). The generalized characteristic can be 
described by the curve -fit equation 


^£ = -0.0435 + 

•2 in 3 

w 10 



1.379 

10 6 



which becomes the following in SI units (Note: In all cases that follow, the second of the 

9 2 

two equations for Ah/Q or for Ap/w is the equation in SI units.); 


= -0. 146 + 


w 


10 " 


86.75 /N 
10 6 



This curve is shown as the solid line in figure 5(b). The dashed line represents a plot of 
points calculated from the equation for the region below the experimental range, from 
N/w = 185 (42. 5) down to N/w = 0. 




Rotational speed, 



Flow rate (of liquid hydrogen), w, Ib/sec 


V 



25 30 35 40 

Flow rate (of liquid hydrogen), w, kg/sec 



(a) Experimental high-speed performance map. 



Speed parameter, N/w, rpm/( I b/sec) 


20 30 40 50 60 

Speed parameter, N/w, rad/kg 

(b) Generalized performance characteristic. 


Figure 5. - Experimental high-speed performance map and generalized performance characteristic for Mark III - 
Model IV turbopump. 
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The pump rotor consists of an inducer followed by radial blades. Pump design data 
and the geometry of the pump are as follows: 


Inlet tip diameter, in. (cm) 6. 95 (17. 65) 

Inlet hub diameter, in. (cm) 1. 65 (4. 19) 

Rotor tip diameter, in. (cm) 12. 25 (31. 11) 

Diffuser inlet diameter, in. (cm) 12.68 (32.21) 

Diffuser discharge diameter, in. (cm) 16.28 (41.35) 

Inducer inlet blade angle, deg 7 

Rotor tip blade angle, deg 90 

Diffuser blade angle, deg 9. 25 

Rotor tip port, in. (cm) 0.485 (1.23) 

Diffuser port, in. (cm) 0.550 (1.40) 

Rotor number of blades (discharge) 48 

Diffuser number of blades 8 

Rotor blade blockage (discharge), percent 29.5 

Discharge pipe diameter, in. (cm) 4.83 (12.27) 

Nominal speed, rpm (rad/sec) 22 500 (2356) 

Nominal flow (liquid hydrogen), lb/sec (kg/sec) 75 (34.2) 


Pressure drop characteristic . - The mean blade inlet and exit velocity triangles for 
the design point operation and for zero slip are shown in figure 6. The evaluation of con- 
stants yields 


y = Vg/Vj = 2.26 

k r = U l /V l = 5 - 75 

k s = V V 2 = 7 - 1 

Evaluation of friction losses gives k f = 0. 466 and k f = 0. 070. Using equation (4a) 

i ? r i? ® 

with rg = rg for a diffuser located directly at the discharge of the impeller blades, the 
generalized negative head characteristic is 

-70. 1 + 0. 119 H + M 2 

Q 2 W 10 4 W 
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I 


^ = -2. 67X10 4 + 12. 23 \ 

Q 2 \ Q / 10 4 \ Q / 

The term of equation (4a) containing k f is -1. l(-0. 042X10 4 ) and the term containing 
kj is -0. 6(-0. 023X10 4 ). During startup, the fluid friction effects are very small. To 
simplify analysis, these two terms will be neglected in all future calculations. 



Figure 6. - Velocity triangles for design point operation and no slip 
assumption for Mark III - Model IV turbopump. (All velocities 
are given in ft/sec (m/sec).) 


The generalized negative head characteristic recalculated for liquid hydrogen (speci- 
fic gravity = 0. 07) and written as the pressure drop characteristic in terms of Ap/w 
and N/w is 


10 


2 AP a -11. 15 + 0. 0826 (-) + °- 655 


w 




10 " 



10 2 -37.35 + 1. 198 

w 2 


41. 18 

07 + 10 4 
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The pump characteristic and the calculated pressure drop characteristic (marked 
X = l) are shown in figure 7. The two characteristics do not approach closely at the 
line Ap = 0 as they should. With the assumption of no mechanical losses, this indicates 
that the diffuser entrance losses at the design point are not negligible or that the incidence 
of the m inimum is different from that assumed by taking X = 1. By trial and error, when 
the correction factor X is assigned a value such that K g = Xk g = 0. 45 k g , the result is 
a curve (shown in fig. 7 and marked X = 0. 45) which nearly meets the extrapolated char- 
acteristic at Ap = 0. 

The new pressure drop characteristic is 
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i 


10 2 = -5. 40 + 0. 048 (?) 

• 2 \w/ 

w V / 


10 2 ^P= -18. 1 + 0.697 /^] 


w 


w j 



By equation (5a), the zero-torque or windmilling point is at N/w = 0, that is, at this 
point fluid fed by a pressurized tank will not initiate rotation. Furthermore, the cal- 
culated characteristic lies entirely (except for the point N/w = 0) in the power dissipa- 
tion zone, and the actual pump operation points will be located above the calculated curve. 
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Figure 7. - Determination of startup characteristic for Mark III - 
Model IV turbopump. 


Comparison of experimental and calculated data. - Two types of tests were conducted 
on the Mark in - Model IV turbopump. In the first, liquid hydrogen was forced through 
the pump using only tank pressure. No power was applied through the turbine. In the 
second type of tests, high pressure hydrogen gas powered the turbine and the turbine dis- 

n 

charge to a subatmospheric vent at 12 pounds per square inch absolute (8. 3 N/cm abs). 

A speed control loop manipulated the turbine power control valve to hold speed constant. 
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For both types of tests, a remotely controlled pump discharge valve was used to load the 
pump. 

The results of the first type of tests are shown in figure 8; no rotation was initiated 
in spite of strong pressurization of the run tank. The curve fit of the experimental data 
for pressure drop across the pump during the zero rotation test is Ap g /w = -0. 057 
(-0. 191) and Ap/w 2 = -0. 545 (-0. 183). 



5L- 


V 

0 


2 4 6 8 10 12 
Flow rate, w, Ib/sec 


0 1 2 3 4 5 

Flow rate, w, kg/sec 

Figure 8. - Pressure drop at zero speed for Mark III - 
Model IV turbopump. 


The results of the second test are shown in figure 9. The setup with the turbine con- 
nected to the exhaust system allowed the leakage from the pump bearing to expand in the 
turbine rows. Some torque was supplied to the impeller to drive the pump, and rotation 
was initiated. The experimental data show a gradual decrease in pressure drop across 
the pump as speed increased. During several runs, tank pressure was varied from 35 
to 50 pounds per square inch (24. 1 to 34. 5 N/cm 2 ) with no apparent influence on startup 
performance. It should be noted that the experimental points in the negative Ap-zone 
cannot be connected by a simple segment of a parabola. 

In the same figure, the calculated curves from figure 7 are shown, and the relation 
between the predicted operation points and the experimental data is displayed. The zero- 
speed point pressure drop parameter is Ap/w = -0. 0545 (-0. 183) while the calculated 
Ap/w 2 = -0. 054 (-0. 181). 
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Speed parameter, N/w, rad/kg 

Figure 9. - Verification of analytically established characteristic for pressurized startup of Mark III - Model IV turbopump. 


Figure 9 is similar to figure 15 (refer to appendix B) where the pressurized startup 
test data and calculated characteristic of an axial pump are shown. In both figures, there 
is a discontinuity of the slope of the pump and pressurized startup characteristics; in 
both cases, the pressure drop and the extrapolated pump characteristics approximately 
meet at the line Ap = 0. Also, the actual zero -speed points are located below the zero- 
speed points obtained by extrapolation of the high-speed characteristic. 


Mark 29F Turbopump 

Data . - The predicted and the experimental high-speed map of the Mark 29F liquid 
hydrogen turbopump (unpublished data from the Rocketdyne Division of North American 
Rockwell Corp. ) is shown in figure 10(a). The generalized characteristic calculated by 

. . ft 

means of the experimental data, in terms of N/w and Ap/w , is 
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This curve is shown as the solid line in figure 10(b) . The dashed line represents the 
analytically extrapolated high-speed characteristic below the experimental range down 
to N/w = 0. 


Predicted performance 

O Experimental data 



till 11 J 1 

0 .1 .2 .3 .4 .5 .6 .7 

Flow rate, Q, nvlsec 


(a) Experimental performance map. 

Figure 10. - Experimental performance map and generalized high- 
speed and calculated pressure drop characteristics for Mark 29F 
turbopump. 


Pump design data and more important geometry data (unpublished data from the 
Rocketdyne Division of North American Rockwell Corp. ) are as follows: 
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Inducer inlet tip diameter, in. (cm) 8. 15 (20. 70) 

Inducer inlet hub diameter, in. (cm) 3. 20 (8. 13) 

Inducer discharge tip diameter, in. (cm) 9.00 (22.86) 

Inducer discharge hub diameter, in. (cm) 7. 50 (19. 05) 

Impeller inlet tip diameter, in. (cm) 6.80 (17.27) 

Impeller inlet hub diameter, in. (cm) 4. 50 (11.43) 

Impeller tip diameter, in. (cm) 11. 50 (29. 21) 

Diffuser base circle diameter, in. (cm) 12. 20 (31. 00) 

Diffuser discharge diameter, in. (cm) 14. 50 (36. 80) 

Inducer mean flow angle at exit, deg 30 


Diffuser mean blade angle at discharge, deg 

Impeller tip blade angle, deg 

Diffuser blade angle at discharge, deg . . . 


Nominal speed, rpm (rad/sec) 28 000 (2932) 

Nominal flow, gpm (m^/sec) 10 HO (0. 638) 



0 20 40 60 80 100 

Speed parameter, N/w, rad/kg 

(b) Experimental generalized high-speed and calculated 
pressure drop characteristics. 

Figure 10. - Concluded. 
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Because the configuration of this pump is very special, its cross section is shown in 
figure 11. The pump consists of an inducer followed by a long diffuser which feeds a 
centrifugal impeller. The impeller discharges into a radial diffuser. The inducing stage 
has a two -row rotor and stator. This special feature is of no consequence in the present 
analysis which considers only the inlet and outlet velocity triangles. 

Pressure drop characteristic. - The velocity triangles for the inducer and the centri- 
fugal stages are established using unpublished data from the Rocketdyne Division of North 
American Rockwell Corporation. The negative head characteristic for the inducer stage 
is calculated using the equation reported in reference 1. The design point deviation is 
assumed, and friction losses are neglected. The result is 

-45. 3 + 0. 0557 fiiV 0,0105 M 2 

q2 \Q/ 10 3 \Q / 

_1.78X10 4 + 5.72M + 2^/nV 

Q 2 W 10 3 W 

The centrifugal stage negative head characteristic is calculated using equation (4a) and 
the assumption of zero slip. The calculated head drop includes the discharge losses. 

The result is 



-82.80 + 0. 1035 


*n\ + 0, 0338 /n\ 2 

&) 10 3 \ Q / 


— = -3. 15X10 4 + 10.65 M /-V 

Q 2 \ Q / 10 3 VV 


The flow in the long diffuser between the inducer and the centrifugal impeller need not be 
considered because only total head changes are calculated and the friction losses are 
neglected. 

Adding the partial results and converting Ah and Q into Ap and w using a speci- 
fic gravity of 0. 07 for liquid hydrogen yield the following for the Mark 29F pump: 
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This curve does not meet the extrapolated high-speed characteristic at Ap = 0. Applying 
the correction factor X = 0. 88 gives the new pressure drop characteristic as 


Ap _ 
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0. 99 /N) 

10 3 W 


+ 1- 34 /W\ 2 

+ 10 6 W 


= -0. 57 + 
w 2 


14.3 

10 3 



84.4 
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This curve does meet the extrapolated high-speed characteristic as required. By equa- 
tion (5a), the zero -torque or windmilling point of the rotor assembly is located at N/w 
= 52. 4 (12. 1). Evaluated separately for the inducer, it is N/w = 99 (22. 9); and for the 
centrifugal impeller alone, it is N/w = 30. 3 (7. 0). These results are in agreement with 
the type of rotor considered. A nearly axial blading will exhibit a high N/w number, 
while a nearly radial centrifugal impeller will have a much smaller zero -torque point. 

Comparison of experimental and cal culated data . - Startup experiments on the Mark 
29F pump were done (unpublished data from the Rocketdyne Division of North American 
Rockwell Corp. ), and the performance data are reproduced in figure 12. The two curves 
give the pressure drop for a range of flows at windmilling and at zero speed (locked 
rotor). Because the tests were run with liquid hydrogen at several temperatures, the 
density of the hydrogen varied, and the pressure drop parameter Ap/w is not a con- 
stant. For the locked rotor tests, it varies between -0. 164 and -0. 171 (-0. 55 and 
-0. 574); and for windmilling test, it varies between -0. 107 and -0. Ill (-0. 359 and 
-0.372). 

Figure 13 gives a comparison of the calculated and the experimental data. The calcu- 
lated curve with the correction factor X = 0. 88 meets the extrapolated high-speed char- 
acteristic on the Ap = 0 line as required. The experimental locked rotor pressure drop 
data verify the calculated point. The windmilling experimental data are located only 
slightly below the corresponding point given by the calculated characteristic and show 
generally good agreement with the calculated curve over the limited range available. 
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Figure 13. - Verification of analytically established characteristic for 
pressurized startup of Mark 29 F turbopump. Test fluid, liquid 
hydrogen. 
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M-l Oxidizer Turbopump 


To illustrate more completely the use of the developed equations and to verify the 
analytical approach, the pressure drop characteristic was calculated for a vaneless 
diffuser pump. The pump is a 3/8 scale model of the M-l liquid oxygen centrifugal pump. 
This scale pump had been tested at high speed in water, but there was no information 
available on its startup performance. Consequently, the conclusions of this study can be 
qualitative only. 

Data. - The M-l scale pump has an axial flow inducer followed by radial vanes and 
its impeller is equipped with back vanes for axial thrust control. The experimental head 
characteristic of this pump in cold water from unpublished data of Aerojet -General Cor- 
poration is shown in figure 14(a). This characteristic corresponds to the following pres- 
sure drop characteristic: 


10 3 -2.478 

w 2 



+ 0. 01565 



10 3 = -8.3 - 0. 169 + 0. 98 


4xl(T 5 



.05 .10 .15 .20 .25 . 30x10 ~ 3 

Flow parameter, Q/N, m 3 /rad 

(a) Experimental generalized characteristic. 

Figure 14. - Experimental generalized and 
analytically established pressurized startup 
characteristics for 3/8 scale model M-l pump. 
Test fluid, water. 
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Total pressure parameter, Ap/w 2 , (N/cm^)/(kg/sec)^ 


This characteristic, shown as the solid line in figure 14(b), has been extrapolated beyond 
its experimental range, down to N/w = 0, as shown by the dashed line. 
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(b) Analytically established characteristic of pressurized startup. 

Figure 14. -Concluded. 
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P um p design data and geometry, from unpublished data of Aerojet -General Corpora- 
tion, are as follows: 


Inlet tip diameter, in. (cm) 7. 20 (18. 3) 

Inlet hub diameter, in. (cm) 3. 24 (8. 23) 

Impeller discharge tip diameter, in. (cm) 10. 69 (27. 15) 

Diffuser discharge diameter, in. (cm) 18.72 (47.5) 

Impeller tip port, in. (cm) 0. 855 (2. 17) 

Diffuser port, in. (cm) 0.908 (2.31) 

Impeller vane blade discharge angle, deg 35 

Discharge pipe diameter, in. (cm) ' 4(10.16) 

Nominal speed, rpm (rad/sec) 7480 (783) 

Nominal flow, gpm (m 3 /sec) 2141(8.10) 


Pressure drop characteristic. - The data listed in the previous section are used to 
calculate the loss coefficients and the ratio y = Vg/V^ for no slip operation at the design 
point. The negative head characteristics are calculated using equation (4a) with k g = 0. 
The effect of the back vanes on the pump characteristic are neglected. After several 
intermediate steps of algebra, the pressure drop characteristic is 
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3 = -4. 38 + 0. 08 f ?) + (51 
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10 2 ^2= -1.47 + 0. 116 
w 2 
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The calculated pressure drop characteristic of this pump is shown in figure 14(b) as 

the dashed line marked k = 0. Also indicated in this figure is the zero -torque point 

s 

calculated using equation (5a). 

Comparison of experimental and calculated data. - Examination of figure 14(b) indi- 
cates that the calculated pressure drop characteristic of the vaneless diffuser pump lies 
considerably lower than the extrapolated pump characteristic. As required, they meet 
at the line Ap = 0. No correction of the calculated performance line is required or 
possible because there are no vanes in the diffuser and the coefficient k g = 0. This qual- 
itative evaluation of the startup performance of a vaneless diffuser pump verifies the 
validity of the assumptions on which the diffuser entrance loss coefficient was calculated. 


28 



I 


CONCLUDING REMARKS 

An analytical method for calculation of startup characteristics of centrifugal pumps 
is presented. The hydrodynamic calculation of the pressure drop characteristic (when 
the pump is acting as a hydraulic turbine) is refined by use of an empirical constant to 
account for losses not predictable by analysis. The correction factor is obtained by mod- 
ifying one parameter in the analytical equations until the calculated curve meets or closely 
approaches the extrapolated experimental pump characteristic at the point where the 
change in total pressure Ap = 0. 

The method is verified by test runs on two types of pumps using liquid hydrogen. In 
both cases, the calculated pressure drop characteristics and the calculated zero -speed 
point, when corrected for diffuser entrance losses by the empirical constant, correlate 
satisfactorily with the experimental data. 

The method has been applied also to a vaneless diffuser pump to illustrate the influ- 
ence of the design parameters on the ability to predict the negative head characteristic. 
Because windmilling test data of the vaneless diffuser pump were not available, this 
evaluation is purely qualitative. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 3, 1968, 

731-11-03-03-22. 
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APPENDIX A 


SYMBOLS 


(Symbol lb is used for both pound force and pound mass. ) 


b 

blade height, ft (m) 

s 

blade spacing, ft (m) 

c 

meridional blade length, ft (m) 

u 

tangential velocity at mean blade 

D 

d h 

diameter, ft (m) 

hydraulic diameter, 2sb/(s + b), 
ft (m) 

V 

diameter, Nr/30, ft/sec (m/sec) 

absolute fluid velocity, ft/sec 
(m/sec) 

f 

Fanning friction coefficient 

2 

V 

relative fluid velocity, ft/ sec 
(m/sec) 

S 

acceleration due to gravity, ft/sec 

o 

w 

work, ft -lb (J) 

g c 

(m/sec) 

dimensional conversion factor, 

w 

flowrate, lb/sec (kg/sec) 


32. 17 (lbm)(ft)/(lbf)(sec 2 ) 

y 

proportionality constant, t^/V i 

Ah 

change in total head, ft (m) 

Z r 

number of impeller vanes 

i 

fluid incidence angle, /3g - a g, deg 

Z s 

number of diffuser vanes 

K 

diffuser loss coefficient, 

a 

absolute fluid flow angle with 

K s 

Ah loss V 3 

^s 

T 

respect to tangential direction, 
deg 

k f,s 

impeller frictional loss coefficient 
diffuser frictional loss coefficient 

a 

relative fluid flow angle with 
respect to tangential direction, 
deg 

k r 

impeller inlet loss coefficient 

(3 

blade angle, deg 

k s 

diffuser inlet loss coefficient 

n 

efficiency 

N 

rotational speed, rpm (rad/sec) 

X 

empirical correction factor for 

Ap 

change in total pressure, 

cy 


stator inlet loss 


psi (N/ cin ) 


viscous loss term, ft-lb/lb 

Q 

volumetric flow rate, ft /sec or 


(m-N/kg) 


O 

gpm (m /sec) 

S2 

flow area, ft 2 (m 2 ) 

r 

mean blade radius from axis of 

Subscripts: 


rotation, ft (m) 

a 

axial component 
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d 

design 

u 

tangential component 

e 

pump exit 

1 

rotor inlet 

i 

pump inlet 

2 

rotor outlet 

loss 

refers to loss 

3 

diffuser inlet 

s 

static 

4 

diffuser outlet 
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APPENDIX B 


PERFORMANCE OF AN AXIAL FLOW, LIQUID 
HYDROGEN PUMP DURING STARTUP 

NASA TM X-1213 (ref. 2) covers an analytical prediction and test verification of the 
performance of an axial flow liquid hydrogen pump during startup. Because of close sim- 
ilarity of startup analysis of the axial flow pump and the centrifugal pumps, the data con- 
tained in this appendix will help to understand the approach used in studies of centrifugal 
pumps. 

The axial flow pump analyzed was the Mark IX liquid hydrogen pump. The pump is 
composed of an axial entrance, mixed flow, axial discharge inducer stage, and six identi- 
cal high-pressure axial flow stages (except that the last rotor discharges directly into a 
vaneless volute). The equation for the generalized negative head characteristic is quite 
similar to equation (3). In calculations, the stator loss coefficient k g was not corrected, 
that is, the stator entrance losses were assumed to be negligible at the design point. 

Equation (5a) was used to calculate the zero -torque or windmilling point. Calcula- 
tions were performed using assumptions of the design point deviation and assumptions of 
no deviation on pump blading. 

Tests were performed in a liquid hydrogen facility equipped with a large overhead 
run tank, a turbopump stand, and a hydrogen disposal system. Tests included three 
pressure levels in the run tank. The windmilling tests were followed by bootstrap tests. 

A number of data points were obtained with pump being driven by a gas turbine in both the 
power dissipation zone and the pump operation zone. 


Test Results 

Figure 15 gives the zero -speed (locked rotor) and zero -torque (steady windmilling) 
calculated performance and corresponding test data. The zero -speed data verify the 
analytically established zero -speed curve (using the design point deviation assumption) . 
The zero -torque or windmilling data fell short of the predicted characteristic. 

Figure 16 shows the startup performance data in liquid hydrogen as well as the cal- 
culated characteristic (the dash -dotted line). The calculation was based on the assump- 
tion of design point deviation. The characteristic extending into the pump operation zone 
is the Mark IX pump characteristic. The solid line is experimental, while the dashed 
line is calculated (extrapolated) using a quadratic curve -fit equation of the experimental 
part of the characteristic. The rotation under windmilling conditions actually stabilized 
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at N/w = 192 (44. 3), while the calculated zero-torque point was N/w = 210 (48. 5). The 
difference should be attributed to mechanical loss torque not considered in the analysis. 
The zero -speed point is well verified. General agreement between the theory and experi- 
ment is good. It will be noticed that, in figure 16, the calculated pressure drop charac- 
teristic has a slope different from the pump characteristic and the two nearly meet at the 
line Ap = 0. 



0 1 2 3 4 5 6 

Flow rate (of liquid hydrogen), w, kg/sec 

Figure 15. - Zero-speed (locked rotor) and zero-torque (steady windmilling) theo- 
retical lines and corresponding test data. 
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Total pressure parameter, Ap/w 2 , (N/cm 2 )/(kg/sec) 2 



0 10 20 30 40 50 60 70 80 90 100 

Speed parameter, N/w, rad/kg 

Figure 16. - Verification of analytically established characteristic for pressurized startup of Mark IX axial pump. Test 
fluid, liquid hydrogen. 
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APPENDIX C 


LOSSES ASSOCIATED WITH VANED DIFFUSERS IN CENTRIFUGAL PUMPS 


Unpublished data supplied by the Aerojet-General Corporation are given in figures 17 
and 18 to illustrate the problem of evaluating losses associated with vaned diffusers in 
centrifugal pumps. 


Test 

Inumber i 


-017 

-014B 

-015B 

-022 

-023 


7.315 


18. 58 


in. 


7. 75 
7. 75 
7.75 
7.95 
7.95 



Number of 
impeller 
vanes, 

Z r 

24 


Number of 

diffuser 

vanes. 

2 

s 



1 



Figure 17. - Vaned diffuser configuration. 


The Aerojet-General Corporation ran a series of tests on vaned diffusers to obtain 
data for off -design prediction of operation of the Mark III - Model IV turbopump. The 
diffuser vanes had essentially log -spiral mean camber lines from the leading edge to the 
vane overlap (throat) station and straight camber lines from this station to the trailing 
edge. 

Loss coefficients were computed from measured average, upstream and downstream, 
static pressures and calculated velocities. The whirl velocity of the impeller fluid exit 
was assumed to be 90 percent of the wheel velocity, and the flow was assumed to follow 
the vane at the diffuser discharge. 

Figure 17 shows the vaned diffusers used to analyze experimentally the off -design 
losses. This figure contains the geometry of the tested rotor -diffuser assemblies. 
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Test number 


O -017 

□ -014B, -015B 
O -022 



Figure 18. - Test results for vaned diffuser loss 
coefficient. 

Figure 18 gives the experimental results. The tests have been run for positive inci- 
dence only. The minimum loss point cannot be correlated with any particular value of 
the ratio i//3g. The minimum loss coefficient varies from 0. 15 to 0. 32, and this varia- 
tion cannot be correlated with the geometrical data of the tested diffusers. To the right 
of the minimum loss point, the losses steeply increase; while to the left, the increase is 
gradual. From these data, it appears that the off -design performance of these pumps 
cannot be predicted. 
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